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DEFINITION OP SYMBOLS 



Wing span - ft 
Lift coefficient 

Rolling moment coefficient 

Yawing moment coefficient 

Side force coefficient 

Roll due to yaw - 

Damping in roll - 

Damping in yaw - 

Yaw due to roll - 

Indicates differentation with respect to 'T* 
time: c i ( ) 

d(-tAr) 

Non-dimensional moment of inertia about 
the x principal axis 

Non-dimensional moment of inertia about 
the z principal axis 

G-ain of pilot 

Radius of gyration - ft 

Potentiometer gain constant i = 1 , 2,3 ... n 




Mach number 



DSFINITIOL 0? SYi'/IBOLS (continued ) 



Roll rate, deg/sec or radians/sec 

Dynamic pressure - psf 

Yaw rate, deg/sec or radians/sec 

Laplace transform variable 

Time - sec 

Steady state velocity - ft per sec 

Trim angle of attack of principal axis, 
radians 

Sideslip angle - deg or radians 

Aileron deflection - radians or degrees 

Damping ratio of the approximate Dutch roll 
mode depicted by the denominator term of the 
"pilot gain" root locus equation 

Damping ratio of short period Dutch roll mode 

A root of the characteristic equation 
(Subscript s refers to spiral root and 
subscript r refers to roll root) 

Airplane density factor - oS b 

Density of air - slugs/ft^ 

-m 



Time parameter pgy - sec 
Roll angle - deg or radians 



Yaw angle - deg or radians 
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DEFINITION OF SYT1BOLS (concluded) 






COrujj 



Undamped natural frequency of approximate 
Dutch roll mode depicted ty the denominator 
term of the “pilot gain" root locus equations 
radians/ sec 

Undamped natural frequency of the short 
period Dutch roll mode - radians/sec 



Subscripts : 

1. The subscripts |3 , and 6a indicate the partial 
derivative of the coefficient with respect to 
the specific subscript. 

2. The subscripts x and z indicate reference axis. 

3. A dot above a variable indicates a derivative 
with respect to time. 

Other terms or symbols are explained where introduced. 
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AESTRACT 

A physical description of the destabilizing 
mechanism responsible for certain pilot induced 
lateral -directional instabilities is included 
herein. In addition, the theoretical and experi- 
mental effort propcrting to support the physical 
arguments is included. 

The basis for this study was a pilot induced 
instability known to exist for the X-15 research 
aircraft-pilot combination for certain conditions 
of flight. Accordingly, a major portion of the 
effort was concentrated upon the analysis of a 
specific known instability. Nevertheless, it is 
the opinion of the authors that the resultant 
conclusions are quite general. 






SECTION I 



INTRODUCTION 

The investigation upon which this report is based 
was concerned with the study of a type of pilot induced 
instability of the lateral -directional short period mode. 

In the more usual case, an instability of this 
nature can be explained as a function of the character- 
istic response time of a human being, in which case the 
problem is quite simply that the human controller can- 
not react to the error signals with sufficient speeds 
Quite obviously, a problem of this nature is closely 
related to the frequency of the mode of motion. 

Some recent studies, however, have indicated the 
possibility that certain pilot induced - or closed 
loop - instabilities can occur which do not appear to 
be strictly a function of frequency. A study by 
Askenas and McRuer, (Ref. 1) based upon root locus 
analyses of a pilot-airplane combination, showed the 
possibility of an instability of this type. This 
study further indicated certain combinations of the 
aerodynamic stability derivatives believed to be 
unfavorable. 

More recently, an analysis was made by Taylor 
(Ref. 2) as a result of simulator studies of the X-15 
research aircraft. These simulator studies predicted 
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certain flight condition boundaries within which the 
X-15 closed loop system would exhibit an unstable 
lateral directional short period mode. (Dutch roll 
mode). That this region did in fact exist was verified 
by flight tests. 

Both of these studies relied heavily on the root 
locus technique for an explanation of how the magnitude 
of pilot initiated control movements could destabilize 
an otherwise stable system. Consequently, although the 
studies showed that the pilot could destabilize the 
system, the actual destabilizing mechanism was not 
explained. This study was initiated in an effort to 
provide this explanation. 

This study was conducted in three separate, yet 
complementary, phases. The theoretical phase consisted 
of an examination of the basic equations of motion and 
the significance of certain parameters upon the resultant 
dynamic behavior or a given pilot-airplane system. Root 
locus analyses formed a major portion of this phase. 

What has been called the computer-simulator phase 
was essentially a study of the same parameters investi- 
gated theoretically, but with the aid of an electronic 
analog computer. In addition, the analog computer was 
used as a fixed base simulator for the investigation of 
an actual closed-loop system. 
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The flight phase consisted of flight tests in- 
corporating the variable stability NAvion aircraft 
described herein. This phase was primarily initiated 
in an effort to provide as realistic a simulation as 
possible of an unstable pilot-airplane combination. 

In general, it may be stated that the results of 
all three phases contributed significantly to the con 
elusions which resulted from this investigation. 
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33CTI0N II 
EQUIPMENT 

The variable stability flight simulator used in 
this investigation was a North American NAvion airplane. 

To provide the variable stability capability, a modified 
Minneapolis -Honeywell autopilot, US\F Type E-12, had 
been installed. Modifications to the autopilot consist- 
ed of the addition of supplementary rate gyros and a 
sideslip transducer. Signals from the transducers, 
proportional to the measured quantities, were introduced 
as input signals to the autopilot, resulting in control 
deflections proportional to the measured quantity. This 
made the system capable of modifying the airplane stabili- 
ty derivatives. Data collection and recording were 
accomplished by means of an A3C0P pulse width, frequency 
modulated telemetering system. An analog computer was 
used in the theoretical development of the problem, and, 
in conjunction with a control stick-operated potenti- 
ometer, provided fixed-base simulation. 

TEST AIRPLANE 

The HAvion is an all metal, low wing, single-engined 
monoplane. The engine, a Continental E-185 , drives a 
variable pitch Hartzel propeller and is rated at 185 
horsepower - maximum continuous at sea level at 2300 rpm. 
The control surfaces are of conventional design. A 
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streamlined static balance is fixed at the outboard end 
of each aileron. The trim tabs for the aileron and 
rudder are of the fixed-bend type, and the elevator 
trim tab is adjustable from the cockpit. The physical 
characteristics of the NAvion are listed in Table I. 

A photograph of the test airplane is shown in Pig. 1. 

The airplane was equipped for this investigation with a 
three-axis variable feedback autopilot. 

AUTOPILOT 

The basic autopilot is capable of altering the 
stability of the airplane in the pitch, roll and yaw 
axes through three channels of operation. For this 
investigation, the pitch channel was rendered inoper- 
ative, with the exception of an unmodified elevator 
control. Each channel is composed of an AC series 
summing network which combines signal voltages from the 
sensing transducers, the electrical flight controls and 
the knob actuated trim controller. The difference be- 
tween this combined signal voltage and a voltage feed- 
back signal from the servo-drum position transducer, an 
error signal, is applied to an amplifier, phase dis- 
criminated and transferred into servo drum rotation by 
means of dual pov/el* relays. The servo drum rotation 
moves the control system of the aircraft. The motion of 
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the control surface continues until error signal voltage 
at the input to the amplifier is reduced to a value 
below the threshold of the system. This value is about 
fifty millivolts. 

The basic autopilot had been modified by deleting 
some of the original provisions and by supplying ad- 
ditional feedback loops. The directional coupler feed- 
back loop, the flight control stick trim provisions, and 
the coordinated turn provisions had been eliminated. 
Sideslip. angle and roll rate feedback loops and a cross- 
control signal proportional to aileron control stick 
deflection had been added to the yaw channel. Sideslip 
angle and yaw rate feedback loops had been added to the 
roll channel. The manually adjustable ratio potentiometers 
were relocated on a gain control panel placed between 
the pilots' seats. An electronic safety device was 
installed to prevent "hard-over" signals by automatically 
disengaging the autopilot if an error voltage exceeding 

r 

approximately 3 volts should instantaneously appear at 
the amplifier input. 

Quantities to be telemetered were taken from the 
appropriate terminals on the main terminal board to a 15 
channel f ilter-conversion unit which converted the voltage 
to DC. A gain control for each channel permitted the 
DC voltage to be adjusted to that required by the tele- 
metering system ( 0-5 volts). Prom the filter-conversion 
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unit , the DC voltages representing the measured quanti- 
ties were taken to an auxiliary terminal strip and 
connected to the telemeter transmitter unit. 

The electrical controls are located in the right 
seat of the NAvion, and the standard manual controls in 
the left seat. Reference 3 described a spring centering 
device and fluid dashpots which had been incorporated 
in the electric stick. These had provided objectionable 
feel characteristics, and , according to the recommendations 
in Reference 3, were removed prior to the flight test 
phase of this investigation. 

INSTRUMENTATION 

The transducers in the feedback loops of the modified 
autopilot provided signals to the autopilot amplifier 
through the AC summing networks and provided the signals 
required to measure the aircraft motions. The flight 
condition data (airspeed, altitude, and outside air 
temperature) were obtained from standard aircraft 
instruments. Errors in these instruments were assumed 
negligible. The physical characteristics of the com- 
ponents of the instrumentation system are listed in 
Table II. 

Two of the modified rate gyros were located on the 
autopilot chassis on the main equipment table and two 
were mounted on the floor of the aircraft in the equip- 
ment compartment o 
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The use of a single sideslip vane was possible 
because it contained two separate potentiometers. The 

i 

sideslip vane was installed on a boom extending four 
feet ahead of the wing leading edge at the wing tip. 

It was assumed that this distance was adequate to mini- 
mize measurement errors due to the wing pressure field. 
The position of the sideslip vane was such that the ya wr- 
ing rate of the aircraft about the vertical body axis 
would affect the sideslip angle measurement slightly, 
but this effect was neglected. 

The telemeter transmitter unit consists of a 
rotary switch sampling each of 43 input channels plus 
two synchronizing channels at the rate of 20 times per 
second. The sampled data are converted to pulse width 
form by a keyer unit and transmitted as a UHF frequency 
modulated signal to the telemeter ground station. The 
filter-conversion unit limited the data measurements to 
14 quantities plus a full scale reference voltage. By 
jumper wiring on the auxiliary telemetering terminal 
strip it was possible to sample a given quantity more 
than once for each revolution of the switch. Specifi- 

I 

cations of the telemeter transmitter unit are listed in 
Table II. 

GROUND STATION 

The ground station data collection equipment con- 
sisted of an A3C0P M series PW/FM ground station, an 
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Arapex Model 309C dual track tape recorder and associated 
graphical recorders. 

The ASCOP ground station received the modulated 
signals from the airborne unit, demodulated and decoded 
the signals to provide for each of the 43 channels of 
information, and provided a continuous voltage output 
to the recorders representing the measured flight data. 

The tape recorder had provisions for simultaneous 
recording of telemeter and voice transmissions. This 
permitted the recording of description of runs, pilot's 
comments and ground station operator's comments along 
with the telemeter data. 

Sanborn four channel pen recorders, Model 154-lOOB, 
were used in the graphical recording phases. 

ANALOG- COMPUTER 

The analog computer used in this investigation was 
a Goodyear Aircraft Corporation Model L3 (GEDA) linear 
electronic differential analyzer. Twenty-four auto- 
matically stabilized DC computing amplifiers were avail- 

7 

able with open-loop gain greater than 5 x 10 , and of 
negligible drift. The computer incorporated an auto- 
matic error indicator and had a guaranteed accuracy of 
one percent. Provisions were available for accurately 
setting computer board potentiometers by the use of a 
special calibration potentiometer and null indicator. 
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Attenpts to "fly" the computer were conducted "by 
means of the simulator shown in Figure 2. This consisted 
of a conventional control stick which, when deflected, 
varied the voltage at the <^ q terminals of the analog 
computer circuit by changing the setting of a potenti- 
ometer. 




Analoq Computer - Simula-hor 
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SECTION III 

PROCEDURE » 

A, In-Flight Simulation 

In-flight simulation of the pilot induced lateral- 
directional oscillation was accomplished with the 
variable stability NAvion. It was expected that results 
obtained in this manner would provide a more suitable 
basis for a physical explanation of the phenomenon than 
results obtained from a fixed base simulator, since the 
full range of visual and kinesthetic motion cues would 
be present. 

Reference 2 points out uncontrollable combinations 
of Mach number and angle of attack encountered with the 
X-15 research aircraft with lateral stability augment- 
ation off and the pilot attempting to control bank angle 
in a normal manner with ailerons alone. Using these 
conditions as a guide, values were obtained for the 
dimensional stability derivatives to be reproduced in 
the NAvion. These are listed in Table III. The calcu- 
lation procedure is shown in Appendix A. These deriva- 

f 

tives were dimensionalized for a dynamic pressure of 
153 lbs. /ft. , corresponding to one g flight. This 
choice of a relatively low dynamic pressure was made in 
order that the capabilities of the autopilot/airplane 
system not be exceeded; it was anticipated that the 
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frequency of the Dutch roll oscillations would be 
reduced, while the character of the pilot-induced 
instability would be retained. Simulation at higher 
dynamic pressure was confined to the fixed base simulator. 

Dimensional derivatives for the NAvion were 
obtained for 6500 ft. standard density altitude and 
120 miles per hour true airspeed (Ref. 4). Autopilot 
feedback gain constants were then determined which would 
alter these derivatives, making them equal to those 
chosen as typical for the X-15. This procedure is shown 
in Appendix A. 

Feedback gain constants of the autopilot could be 
altered by adjustment of potentiometer settings in the 
cockpit of the NAvion. In order to determine gain 
potentiometer settings which would provide the required 
constants it was necessary to calibrate the system.' 

Bach autopilot rate gyro was calibrated by removing 
it from the airplane and remounting it on an electrically 
driven turntable which could be rotated at a controllable 
rate. In this manner control surface deflections pro- 
duced by given angular velocities could be established. 

The deflections were measured by a propeller protractor 
for various gain potentiometer settings. 

The sideslip vane was calibrated by moving it through 
a specified angle and simultaneously measuring the control 
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surface deflections produced at different gain potenti- 
ometer settings. Electric stick calibration was obtain- 
ed in the same manner, a propeller protractor being 
used to measure stick deflection. 

The results of the calibration are shown in Figure 
3, as curves of Feedback Gain Constant versus Gain 
Potentiometer Faceplate Setting. The numbering of the 
constants corresponds to the number of the potentiometer 
by means of which each may be adjusted. The cockpit 
panel arrangement is shown in Figure 4. 

In order that a quantitative interpretation of the 
telemetered motion histories of the aircraft could be 
made, it was also necessary to calibrate the system 
for telemetry. This was performed in exactly the same 
manner as the calibration described above, with the 
exception that the measured output quantity in this 
case was voltage at the telemetry transmitter package, 
rather than control surface deflection. It was decided 
that the following quantities would provide sufficient 
information for the analysis to be performed: roll rate, 

yaw rate, stick deflection and sideslip angle. The 
results of the telemetry calibration are shown in 
Figure 5 as curves of tne foregoing quantities versus 
Telemetry Input Voltage. 
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Longitudinal feedbacks, (pitch angle, pitch rate, 
velocity and angle of attack) were not employed in this 
investigation. The elevator gain potentiometer was 
adjusted to a value which gave a satisfactory elevator 
response. 

Certain deficiencies existed in the equipment which 
made it difficult to achieve an exact simulation. Among 
these were control cable stretch, structural elasticity 
of the aircraft, inaccurate resolution due to finite 
autopilot "peck sizes," very small potentiometer settings 
and the difficulty of accurately reproducing a potenti- 
ometer setting. As a result of these deficiencies it 
was found that the simulation procedure outlined in this 
section provided only a first approximation to the desired 
aircraft response. It was necessary, therefore, to 
further adjust the cockpit gain potentiometer settings 
in order to achieve the desired results. 

This arbitrary correction to the calculated potenti- 
ometer settings would necessitate matching of aircraft 
responses to analog computer responses if the objective 
was exact simulation. A reasonably good qualitative 
simulation of such motions can be readily achieved, 
however, and was considered to be sufficient for the 
purpose of this investigation. 
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E. Analog Computer Study 

Three basic tools were used to accomplish the 
computer phase of this investigation. These were: 

1. the G^vDA Analog Computer , 

2. the Sanborn Recorder, and 

3. the Control Stick Simulator. 

The computer and recorder were used to do the open loop 
studies of the system, while all three units were used 
to complete the closed loop work. 

The equations mechanized in the computer are listed 
and discussed in detail later in this report. 

When open loop response to step inputs by the 
computer matched calculated values as to time to damp 
to 1/2 amplitude and period of oscillation, attempts 
were made to "fly" the computer. These closed loop 
trials were accomplished by using the control stick 
simulator in the same manner as the control stick in an 
aircraft. The time histories of bank angle, roll rate, 
sideslip angle, and yaw rate on the recorder were used 
as references. The control stick simulator produced 
input voltages to the computer analogous to aileron 
inputs to the basic airplane. 

Three separate configurations of the basic equations 
of motion were investigated. These were: 
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A. 0d o = 0 and low dynamic pressure (153 psf). 

B. Od 0 = 10° and low dynamic pressure (153 psf). 

C. <X 0 = 10° and high dynamic pressure (1000 psf) 

The wiring diagrams for the computer study are 



shown in Figures 6 and 7. 




FIG. (o 

ANALOG COMPUTER DIAGRAM 
Lateral - directional Equations q = 153 lb. /ft. 




FIG. 7 

ANALOG COMPUTER DIAGRAM 

Lateral - directional Equations q = IOOO lb. / ft. 2. 
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S3CTI0K IV 

DISCUSSION and results 
T heoretical Development : 

During the past several years, considerable emphasis 
has been placed on the determination of what parameters 
can best be used to measure the lateral-directional 
handling qualities of an airplane. As a result, a consider- 
able body of literature has come into being and certain 
airplane response characteristics have come to be re- 
garded as unfavorable or undesirable from the standpoint 
of control by a human pilot. (Ref. 1). Primarily as a 
result of the human pilot being involved as a part of the 
feedback control loop, investigations of this nature have 
not been easy to interpret either physically or mathe- 
matically. 

This particular investigation has been concerned 
with a rather special segment of the lateral-directional 
handling qualities problem - specifically the case where 
the inclusion of a pilot in the control loop has resulted 
in an unstable Dutch roll oscillation, although the open 
loop Dutch roll response of the aircraft was known to be 
stable. In short, this investigation has been concerned 
with the specific problem of a pilot induced Dutch roll 
instability. 

The problem of a pilot-airplane combination of the 
nature described above could have been regarded as general 
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and an overall solution sought. However, it was felt 
by the authors that a more concentrated effort could be 
made if the problem were approached from the opposite 
point of view - namely a detailed theoretical and 
experimental investigation using, as a guide, a pilot- 
airplane combination for which the results were known. 

The reasoning followed was that a detailed study of a 
particular situation might then yield a theory which 
could be investigated in a more general sense. 

According to Ref. 2, the X-15 research airplane- 
pilot combination was found to result in an unstable 
lateral-directional oscillation for certain flight 
conditions if the stability augmentation system was made 
inoperative. For this reason, the X-15 airplane-pilot 
combination was the closed loop system chosen for the 
detailed study. Specifically, the instability associated 
with a flight Nach number of 3 and a ten degree angle of 
attack was investigated. 

The fact that the above configuration v/as the one 
chosen for study placed several important restrictions 
upon the problem and thereby upon the resultant conclusions. 
First, although the above configuration was handled 
throughout this study as though the aircraft had been in 
steady flight, in reality the X-15 encountered the above 
stated flight condition only in transitory flight. 
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Accordingly, the aerodynamic stability derivatives were 
continuously changing, with the result being that the 
pilots of the X-15 were advised to refrain from using 
the rudder as a control. Therefore, this entire investi- 
gation was conducted assuming the ailerons to be the only 
type of lateral-directional control available , although 
the effect of using rudder is discussed in the flight 
phase section of this discussion. A second condition 
implied by the problem as stated is that the frequency of 
the oscillation was not a factor. That this must be so 
can be reasoned by noting that the frequency is more a 
function of the dynamic pressure than of Mach number, 
while the instability has been described as being a 
function of Mach number and angle of attack. Thus, the 
solution sought was not one dependent upon the inability 
of a human being to respond with sufficient speed to an 
error signal. In effect, a suitable explanation to the 
problem encountered by the X-15 would have to be valid 
even if the pilot were assumed to respond without a time 
lag. 

The stability derivatives required for a detailed 
study of the above condition were obtained from Ref. 5 
and are listed in Table IV. These derivatives lead to 
autopilot potentiometer gain constants listed in Table V. 

It is noteworthy that only two of the stability derivatives 
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were significantly different from what one might ordi- 
narily expect. Namely, hothC n ^ a , (aileron yaw) and Cfi.jo, 
(dihedral effect) had a positive sign. Thus, these two 
derivatives immediately became suspect and it was expected 
that a theoretical analysis of the problem would show 
either one or both of these derivatives to be involved 
in an explanation of the instability. In this rather 
special case of a single control, pilot-induced insta- 
bility, the unusual sign of the aileron yaw became of 
particular interest. 

A proposed theoretical explanation to the problem 
encountered by the X-15 aircraft was published by Taylor 
(Ref. 2). In his work, Taylor showed, by root locus 
techniques, that the Dutch roll mode of an aircraft with 
the characteristics exhibited by the X-15 could be made 
unstable if a pilot gain proportional to roll angle and 
roll rate was included in the feedback loop. According 
to Taylor, the ratio of could be used to define 

the region where an instability might occur - a ratio 
greater than one being indicative of an unfavorable 
condition. The human transfer function used by Taylor 
in his analysis was determined empirically and found to 
be K(1 + .57s). Thus a pilot response proportional to 
both roll angle and roll rate was assumed. The analysis 
followed by Taylor was quite believable for two reasons. 
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First, the transfer function was seen to be independent of 
a pilot time lag which would seem to isolate any effect of 
frequency from the basic control problem. Secondly, the 
analysis as carried forth by Taylor resulted in the con- 
clusion that the basic cause of the problem was the positive 
sense of the dihedral effect ( Cj^). What Taylor’s analysis 
did not do was present a physical explanation of the phe- 
nomenon. 

In an effort to assure that the problem had been 
interpreted correctly by the authors, the initial theo- 
retical effort was directed toward matching the results 
of Taylor’s investigation. Thus the human transfer 
function proposed by Taylor was used and a root locus 
study of the effect of varying pilot gain undertaken. 



The set of equations used was as follows: 
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It is worthy of note that the above set is for a prin- 
cipal axes coordinate system. This particular coordinate 
system was chosen for this study inasmuch as the stability 
derivatives as obtained from Ref. 5 had been calculated 
with respect to this system of axes. 

The result of the root locus study can best be 
explained by considering separately the effect of varying 



- 30 - 



the pilot gain, K, and the trim angle of attack, 0<f 0 . 

The fact that 0C o was a significant parameter actually 
was determined during the analog computer phase of this 
program and is discussed in detail in connection with 
that part of the program. Here the fact that <X 0 is import- 
ant will merely he stated and then substantiated by the 
root locus technique. 

Assuming the pilot gain to be the variable of interest, 
the determinant of the coefficients of (1) was expanded 
with the following root locus equation resulting: 

(2) A ~ \ S ) ( \ C + A + U>ny j Aq/, p, 

( \ Z + A + -57 A) 

The derivation of the above equation is included in 
Appendix B. 

Figure 8a was the type of root locus which resulted 
from the above equation when CX 0 was assumed to be negligible,. 
This root locus showed the inclusion of a proportional 
pilot gain to be a "stabilizing factor at zero angle of 
attack, provided the stability derivatives could be 
assumed constant. Furthermore, a study of the parameters 
which determine the value of CO n ^» (Bq. 5 in Appendix B), 
reveals that an increase in either of the two suspect 
stability derivatives (C n & a and C.^) would actually de- 
crease the value of COnp and thus tend to stabilize the 
closed loop system. This is a result of the signs of 
both derivatives being reversed from what might be termed 
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The next step was to include the correct value of C 
(10 degrees). From equations (1) it can be seen that the 
inclusion of OC 0 has no effect upon the denominator term 
of the root locus equation but does alter the character- 
istic equation somewhat. The added terms which are the 
result of including CX 0 are: 



(3) CX c 



£j-j (C I/a C„ r - Cn (3 C{ r )X- -f^X 



r 



J, 



It can be seen from equation (3) that with a positive 
Cg^and all other derivatives of the usual sign, the 
effect of including (X 0 is to decrease the A and A 
terms of the characteristic equation. This had little 



effect upon the roll and spiral modes but did result in 
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both the damping ratio ) and the undamped natural 

frequency (OOmj/ ) of the Dutch roll mode being appreci- 
ably decreased. The root locus which resulted is depicted 
by Pig . 8b . 




Fig 8b. Roo+ Locus "for O C Q - 10° t CKHX 00 

This is similar to what Taylor obtained. Furthermore, 
it can be seen from equation (3) that the destabilizing 
effect of (X 0 can be attributed to the positive sign of 
the dihedral effect. (Cfi^). 

Having shown the importance of the trim angle of 
attack, it was deduced that a root locus using trim 
angle of attack ( OC 0 ) as the varying quantity might be 
informative. This study was carried out by assuming the 
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stability derivatives to be independent of the trim 
angle of attack. Thus only the orientation of the air- 
craft with respect to its flight path was assumed to be 
variable. The root locus equation for this investigation 
had the following form: 

(4 ) (X + \ri(X + \s)( X^ + t is i ~ 

X(X+ K'/K') * a ‘ = 



with the denominator terms found from equation (3) and with 

K. . -/^ an, K- '^(C^Cn, - Cn p C ir ) 

Figure 8c shows the root locus f or 0-O^ o ^00 f or no 
pilot gain while figure 8d shows the same result for an 
arbitrary pilot gain. Calculations for the conditions 
investigated showed the system with pilot gain to go 
unstable for a trim angle of attack on the order of 
twenty five percent of what was required for instability 
without the pilot in the loop. Conversely it was noted 
that for a trim angle of attack of ten degrees, the open 
loop system was so slightly stable that a small pilot 
gain could be expected to be sufficient to cause the 
closed loop instability. From equation (4) it can be 
seen that a negative Cj^ would change the sign of K' 
and in turn change the angle condition for the root locus 



diagram. For this case, the effect of increasing OC 0 
would not be destabilizing. 
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Several interesting observations resulted from a 
study of the root loci briefly discussed above. First, 
it was observed that if the signs of both the dihedral 
effect ( Cfip ) and the aileron yaw (On 6a ) had been’ 
reversed, there would have been virtually no change in 
the closed loop response at low angles of attack. How- 
ever, a study of the cofactor of the ( CX 0 ) term showed 
that a beneficial effect would have been realized as 
the trim angle was increased. Secondly, the observation 
was made that for a positive dihedral effect, the pro- 
verse sense of the aileron yaw definitely tended to 
improve the closed loop stability while an adverse 
aileron yaw would have had the opposite effect. That 
is, an adverse aileron yaw would have placed the complex 
zero (Fig. 8b) even further above the complex pole. 

This was believed to be particularly significant inas- 
much as with pilot time lag removed as a significant 
contributor to the instability, one might reason that 
the only way that the pilot could be contributory as a 
destabilizing factor would be as a means for introducing 
an unfavorable forcing function while actually responding 
correctly to the error signal. In the case of a single 
control, this meant an unfavorable input into either the 
rolling or yawing moment equations by means of an aileron 
deflection. For the particular case investigated, the 
root locus plots showed a proverse aileron yaw (Cn&a ) 
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to be favorable while an increase in the trim angle (CX 0 ) 
was seen to have a destabilizing effect. Thus it was 
reasoned that the probable physical effect was that an 
increase in trim angle of attack opposed the proverse 
aileron yaw effect until, at some angle of attack, the 
sense of the sideslip resultant from an aileron deflection 
was probably reversed. This might be explained by noting 
that a simple uncoupled rolling motion about the body 
axis of an aircraft will result, as the roll angle increases, 
in the conversion of angle of attack into sideslip angle. 

For example, a positive roll would tend to produce a 
positive sideslip. This is the opposite of the effect 
produced by a proverse aileron yaw. Thus it is reason- 
able to assume that at some point the effect of angle of 
attack and of proverse aileron yaw might be equal. At 
this point and beyond, an increase in the pilot gain 
would cause the pilot to introduce a destabilizing force. 
Looked at in the light of Taylor* s analysis, this point 
where the aileron yawing moment due to aileron deflection 
reversed direction is probably coincident with the point 
where becomes greater than unity. 

In summary, the theoretical effect agreed quite well 
with Taylor’s analysis although it is felt by the authors 
that the approach involving the variation of trim angle 
of attack provided a basis for the physical reasoning 
explained above. However, it can be seen that even if 
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the physical effect of increasing the trim angle is as 
stated, the mechanism by which the instability occurs 
cannot be explained from the theory alone. Fortunately, 
a study of the time histories of both the analog computer 
simulator and the variable stability NAvion did provide 
an insight into the mechanism involved and did tend to 
agree with the hypothesis suggested by the theory. 

Analog Computer - Simulator Phase : 

As mentioned previously, the approach used in this 
investigation was to make a detailed study of an airplane- 
pilot combination which results had shown to be unstable. 
For the computer phase of the investigation two variations, 
along with the basic condition, were examined. The basic 
condition v/as the X-15 research airplane flying at 
M = 3.0 with an angle of attack of ten degrees. The two 
variations were: to keep all other conditions the same 

but decrease the trim angle of attack to zero, and to 
keep all other conditions the same while increasing the 
dynamic pressure from 153 to 1,000 pounds per square foot. 
All three of these conditions wil^L be discussed below. 

Condition I: M = 3.6. q = 153#/ft 2 . CXo = 0 

For this condition, the basic equations of motion, 
using the non-dimensional X-15 coefficients, reduce to 
the following set of computer equations: 
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100(3 = -3.35(3 - 100 0 + 1.0670 
.20 = l .67(3+ .03454- -.0420 + 3.66<$ q 
100 = 60(5 - .451 0 -.0451 0 + 9.96c, 

Both step aileron inputs and "attempts-to-f ly" the 
computer simulator showed this condition to he completely 
manageable without* the slightest tendency to go unstable. 

This fact, in conjunction with the other studies 
outlined herein, confirmed the significance of the trim 
angle of attack ( 0( 0 ) to the study of the closed loop 
instability. 

Condition II: M = 3.0. o = 153#/ft 2 . 04o = 10 degrees 

This condition resulted from an analysis of the 
results of Condition I. By re-evaluation of the basic 
equations, using 0 £ 0 = 10 degrees, the following set of 
equations was obtained: 

100/3 = -3.35(3 - 1000 * 1.0670 + 17-45 0 

10 0 = 90(3 + 1.725 0 - 2.12 0 4- 83.9<S a 

20 0 ^ 55(3 - .9020 -.09020 + 9.08<$ a 

These equations were mechanized into the computer and 

resulted in the predicted pilot-induced instability. 



A step aileron input and an "attempt-to-f ly" the 
computer are shown in Figures 9 and 10. From the step 
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input response there were several important things noted. 
First, it was seen that the motion was lightly damped. 
Second, there was considerable excitation of the Dutch 
roll mode from a coupling of the rolling and yawing modes 
of motion. Third, although the sign of Cn^g was posi- 
tive - leading one to believe the yaw, and consequently 
the sideslip, due to aileron should have been proverse - 
the total effect of the parameters affecting the side- 
slip was just the opposite. That is, the sideslip angle 
response described by these equations, due to a step 
aileron input, was as might be expected from an adverse 
yawing moment. The fourth and last item to be noted was 
the phase angle difference between the angle of roll and 
the angle of sideslip. It could be seen from the time 
histories, that the two motions were almost exactly out 
of phase. The significance of the phasing of these angles 
while the Dutch roll motion was in progress, as well as 
the sign of the sideslip resulting from an aileron input, 
will be discussed more fully later. 

Using the control stick simulator, "attempts-to-f ly" 
the simulator resulted in oscillations of increasing 
amplitude if the ''pilot" responded in a normal manner 
to bank angle and roll rate. Using the stick to make 
aileron inputs in an effort to keep the wings level 
resulted in an increase of the sideslip angle , due to 
the aileron deflection. Thus, sideslip grew, fed the 
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roll excursions, and the resultant situation was one of 
a closed loop instability. By releasing the stick, the 
oscillations damped out, although the roll angle grew 
indefinitely as a result of the unstable spiral mode. 

An "attempt-to-fly" the computer by using the 

technique described by Petersen in Reference 6 was 

completely successful. The oscillations were allowed to 

build up using normal pilot techniques, and then, at 

some point, ailerons were used primarily to stop the 

r 

sideslipping, and secondarily to keep the wings level. 

As the p trace came through zero, the stick controller 
would give the stick a sharp pulse in the direction the 
trace was moving. In the airplane this would correspond 
to a sharp right aileron pulse as the nose of the airplane 
went through zero degrees sideslip heading toward the 
right. In this manner, the motion could be damped out 
within a few cycles. It did not appear feasible to 
maintain steady state conditions using this technique, 
but oscillations resulting from normal pilot inputs could 
be rapidly damped out in this manner. Roll angle control 
was maintained by varying the magnitude of the pulse and 
the stick position betv/een pulses. 

During the study of this condition, the following 
observations were noted: 

1. The oscillatory motion was lightly damped, 

2. There was considerable coupling between 
the roll and sideslip motions. 
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3. The sideslip generated by an aileron deflection 
was adverse even though Cn^a was positive. 

4. During Dutch roll oscillations , the angles of 
bank and sideslip remained approximately 180 
degrees out of phase. 

5. Normal atternpts-to-f ly were unsuccessful, and 

6. The j5 technique could be used to damp out any 
Dutch roll oscillations that developed. 

Condition III: M = 3.0* 0 = 1000#/ft^. Cy 0 = 10 degrees 

This condition was investigated primarily to see if 
the characteristics of the motion were affected by changes 
in dynamic pressure. The basic equations of motion modi- 
fied to increase the dynamic pressure were written as: 

20(3 = -4-52(3 - 200 + 1.440 -t- 3-490 
20 = I 16.6(3 +2.32 0 - 2.840 +!08.4<Sq 
40 = 7i .2(3-1.2120 -. 12120 +K.72&Q 

These equations were mechanized in the computer and 
both step aileron inputs and "attempts-to-f ly" were 
observed. The basic character of the motion did not 
change. All the comments recorded under Condition II 
are valid for the higher dynamic pressure., The only 
significant change in the motion was the frequency of 
oscillation of the Dutch roll mode. Previously the 
motion had a period of about five seconds. At this in- 
creased dynamic pressure the period was reduced to 
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approximately two seconds. This did not alter the diffi- 
culty of attempting to fly the configuration. Nor did 
the increased frequency invalidate the fb technique of 
control. 

This condition can be summed up simply by saying it 
was identical to Condition II except things happened a 
little faster. However, the motion was not so fast as 
to make the controller’s task too difficult, nor so fast 
that the character of the motion and subsequent responses 
to the motion were changed. 

One result of the above observations was that the 
frequency of the motion (v/ithin the limits investigated) 
was seen to be an insignificant contributor to the in- 
stability. As long as the motion proceeded at a pace 
which the controller could follow, and as long as the 
normal control motions were made, the basic instability 
occurred. Thus the frequency of the motion was elimi- 
nated as an important parameter insofar as a study of 
the basic destabilizing mechanism was concerned. This 
confirmed the hypothesis outlined in the theoretical 
portion of the discussion. 

Flight Phase : 

The flight phase of this investigation involved the 
use of the variable stability NAvion. The physical 
characteristics of this airplane, its simulation equip- 
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men t and the accompanying instrumentation are described 
elsewhere in this report. 

The initial flight phase was based upon an in- 
stability encountered by the X-15 research aircraft in 
a steady flight (lift = weight) condition at Mach 3 and 
at ten degrees angle of attack. The dimensional X-15 
derivatives used in calculating the potentiometer 
settings are listed in Table III. A direct method for 
simulating the X-15 side force equation was not available 
and the dynamics introduced by the autopilot equipment 
were not known by the authors, hence it was necessary to 
correct the computed autopilot potentiometer gain 
settings in flight until a desired condition was achieved 

Flight simulation was directed toward achieving the 
desired character of pilot induced instability, rather 
than toward an exact simulation of the dynamic response 
of the X-15. On this premise, and on the basis of 
computer studies, it was decided that the most signifi- 
cant characteristics to be simulated were the following: 
1) a lightly damped Dutch roll mode 2) a roll rate 
response to an aileron function which contained a signifi 
cant Dutch roll effect. 

Cockpit potentiometer settings were adjusted until 
the two above characteristics were obtained. It was 
found that the period of the Dutch roll mode obtained 
in this manner was approximately three seconds, and the 
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steady state roll rate of the aircraft was slightly de- 
creased. Thus Dutch roll period is of the same order 
of magnitude as the periods of the Dutch roll modes 
studied on the computer. In addition, computer studies 
indicated that the frequency of the Dutch roll mode was 
unimportant insofar as the fundamental problem was 
concerned . 

Undoubtedly the most important characteristic was 

I 

simulated, namely, the closed-loop pilot-airplane system 
exhibited an unstable Dutch roll oscillation in spite of 
the fact that the open loop response was lightly damped. 
Thus, as in the case of the X-15, the attempts of a 
pilot to improve upon a marginal condition had a de- 
stabilizing effect. Significantly, this unstable result 
could be avoided if the pilot used both ailerons and 
rudder for control. Just why the authors feel this 
point to be significant will be discussed later in this 
report . 

Time histories of the various angles and rates as 
recorded on a typical flight are included herein. In- 
asmuch as the mode of interest was the Dutch roll , the 
authors felt the two most significant time histories to 
be those of the sideslip angle (j3) and the roll rate 
( 0 ). Accordingly, the time histories of these two 
items were studied in detail and will be the items 



discussed herein 
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It can 'be seen from the time histories of the air- 
plane response to a step aileron input (Fig. 11) , that 
the roll rate did contain a considerable amount of Dutch 
roll, and further, that the Dutch roll motion was lightly 
damped. Also, it can be seen that the proverse sense of 
the aileron yaw resulted in the initial sideslip generated 
by an aileron input being of opposite sign to the roll 
rate response. This was as might be expected for an air- 
plane with appreciable proverse aileron yaw and only a 
slight amount of yaw due to roll. However*, this initial 
sideslip response to an aileron input was the opposite 
of what was observed to happen on the analog computer 
simulator. With the benefit of hindsight, it can be 
stated that this difference was probably due to the fact 
that the effect of the trim angle of attack (Cy 0 ) had 
not been considered. 

From the time histories of the attempts to fly the 
resultant configuration, (Fig. 12) it can be seen that 
attempts by the pilot to control the oscillation resulted 
in excursions of increasing magnitude for both the side- 
slip angle and the roll rate. Releasing the controls, 
(Fig. 13) on the other hand, resulted in decreasing 
amplitudes of motion. However, once large amplitudes of 
motion had been obtained, what appeared to be a limit 
cycling of the autopilot prevented a complete damping of 
the open loop motion. An important piece of information 
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that can be obtained from these time histories is the 
phase relationship of the sideslip angle ((d) and the 
roll angle ( 0 •) - the latter being assumed to lag the 
roll rate by a quarter of a cycle. In particular, it 
can be seen that the sideslip angle and the roll angle 
are almost exactly in phase. This was also in contrast 
to the phasing obtained by the electronic analog simulator. 
Thus it was noted that both the analog simulator and the 
NAvion simulator exhibited a pilot induced oscillation 
in spite of the fact that the character of the motion was 
virtually opposite for the two cases. That the sideslip 
resulting from a step aileron input was also different 
for the two systems suggested the possibility that, in 
spite of the apparent difference, the same destabilizing 
mechanism was involved in each case. Accordingly, the 
aileron yaw was reversed from the proverse to the adverse 
sense and the tests repeated. Time histories (Pig. 14) 
show that although the response of sideslip angle to an 
aileron input was reversed, the percent of Dutch roll 
in roll was maintained at a high level and the Dutch 
roll remained lightly damped. However, this configuration 
could be controlled by the pilot. This is not to say 
that the resultant configuration was a good - or even 
acceptable - airplane from a handling qualities stand- 
point. The damping of the Dutch roll mode was so slight 
that unpleasant oscillations developed rather easily - 
particularly if rapid control reversals were attempted. 
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However, it is significant that the pilot was able to 
control the aircraft, and, in fact, damp the oscillations 
when they occurred using only normal aileron control 
movements. This seemed to confirm the hypothesis that 
the destabilizing mechanism for the X-15 pilot-aircraft 
combination (as manifested by the analog simulator) was 
of a similar character to what was obtained with the 
NAvion. In particular, it could be seen that the in- 
stability resulted when the pilot, in an attempt to 
control roll angle, succeeded in reducing the roll angle 
only at the expense of increasing the sideslip angle. 

During this flight phase, flight of the basic 
(closed loop unstable) configuration using both ailerons 
and the rudder to maintain lateral-directional control 
was tried. The results showed that by using the rudder, 
the sideslip angle could be controlled and the aircraft 
handled comparatively easily. Rapid reversals presented 
a problem as it was difficult to keep the aircraft in 
balanced flight throughout the maneuver and a rather 
large Dutch roll oscillation often resulted. However, 
by using rudders, the oscillations could be damped within 
a few cycles, (Pig. 15). This also tended to confirm the 
importance of the sideslip angle to the oscillation. 

Aerodynamic derivatives other than the aileron yaw 
were varied and a pilot opinion of these variations 
obtained. However, even though many of these derivatives 
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v/ere found to have an important effect upon the handling 
qualities, none were found to he so directly involved in 
the actual destabilizing action of the pilot as the 
aileron yaw, but rather were indirect contributions to 
the problem. As such, the effects of these other stability 
derivatives were probably of importance in establishing 
the conditions believed to be necessary for a problem 
of this type to exist. These conditions are further 
discussed in the following section. 

The reason for the difference in phasing of roll 
and sideslip between the computer and h Avion is believed 
to lie in the inability to directly simulate the side 
force equation with the hAvion. An attempt was made to 
adjust the NAvion stability derivatives in such a way as 
to allow for the effect of trim angle of attack ( 0l o ) . 
Plight test of this altered configuration, however, 
indicated that the phasing had not been changed. 

The fact that the KAvion did not yield the same 
phasing of sideslip and roll as was obtained with the 
analog computer was not considered by the authors to be 
a shortcoming of the flight simulation. It was felt, 
rather, that this difference in phasing helped to 
emphasize the more general nature of the physical 
explanation of the instability mechanism proposed here- 
in - that is, since the same character of pilot induced 
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instability was manifested in both of the two different 
types of systems investigated, the problem is not limited 
to one particular aircraft. Furthermore, since the 
analysis of the mechanism of instability is consistent 
with both of the two types of motion phasing, the appli- 
cations of the analysis also become less limited. 

Summary of Results : 

Quite obviously the different phases of this program 
have resulted in a multi-sided approach to the basic 
problem of a pilot-induced Dutch roll instability. In 
general, it can be stated that each phase of the investi- 
gation contributed something unique. Thus, in this section 
of the report, these separate and distinct contributions 
will be summarized in an effort to show how each con- 
tributed toward the formulation of what is believed to 
be a correct physical explanation of the destabilizing 
phenomenon. 

From the root locus investigations discussed pre- 
viously, it was determined that, with a positive dihedral 
effect, an increase in the trim angle of attack ( (V 0 ) 
had a destabilizing effect upon the characteristic open 
loop response of the X-15. Further, a study of the 
equations of motion showed that the effect of increasing 
the trim angle of attack was to increase the coupling 
of the sideslipping and rolling motions of the aircrafto 



- 57 - 

In addition, a study of the co-factor of the trim angle 
of attack (0£ o ) term, in conjunction with the root locus 
shown in Figure 8c, indicated that the phasing of the 
motions was of such a nature that destabilizing moments 
were generated. Even more significantly, the root locus 
analysis showed that, under certain conditions, the normal 
inputs of a pilot attempting to keep the wings level 
had a destabilizing effect upon the Dutch roll mode. 

It was then reasoned that inasmuch as pilot response 
time had not been considered in the analysis, the de- 
stabilizing effect of the human controller could not 
have been the result of erroneously applied corrections 
to unwanted deviations in roll ai gle , but rather had to 
be due to the indirect effects of the resultant aileron 
deflections. This indirect effect was seen to be the 
conversion of angle of attack into sideslip. Thus it 
was concluded that the attempt of a pilot to control 
roll angle resulted in the generation of a sideslip 
angle. Further, the phasing between the rolling and 
sideslipping motions was of such a nature that one 
motion tended to reinforce the other. 

The results obtained as a consequence of the analog 
computer-simulator phase of this program complemented 
the theoretical investigation and, in addition, provided 
an explanation of the actual destabilizing mechanism. 

The computer studies confirmed the importance of the trim 
angle of attack to the oscillatory behavior of the system. 
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Sub sequent ly the root locus investigation was under- 
taken and the results discussed above obtained. 

Study of the time histories obtained during the 
computer phase of the investigation showed that the 
Dutch roll oscillation, for the equations represent- 
ative of the X-15, was characterized by angles of roll 
and sideslip which were of opposite sign. Further, 
the open loop response to a step aileron input showed 
that an input producing a positive roll resulted in an 
initial positive sideslip (Fig. 9). This was not what 
was expected, inasmuch as the aileron yaw was proverse 
and the yaw due to roll was small. However, as has been 
previously explained, the conversion of angle of attack 
into sideslip as the result of a rolling motion could 
produce an effect of this type. 

By combining the two effects noted above - which 
would be the natural result of a pilot attempting to 
control a Dutch roll motion with aile rons - it could be 
seen that an attempt to control the angle of bank in a 
normal manner would result in an input tending to in- 
crease the angle of sideslip. With the high degree of 
coupling present between sideslip and roll, as evidenced 
by the considerable Dutch roll contribution to the roll 
rate time history, an unstable behavior resulted. This 
hypothesis was checked by artificially increasing tbs 
proverse aileron yaw to such an extent that the initial 
sideslip generated by a step aileron input was changed. 



- 59 - 



That is, an input which resulted in a positive roll 
yielded a negative sideslip. Figures 16 and 17 show 
the time histories resultant from this configuration. 

With the increased proverse aileron yaw, the configuration 
could be controlled in spite of the fact that the Dutch 
roll mode remained lightly damped and the percentage of 
Dutch roll in roll remained appreciable. The lightly 
damped Dutch roll mode did result in easily induced 
oscillations but, unlike the basic configuration, pilot 
attempts to control roll angle did not destabilize the 
system. In fact, as might be expected, it was found that 
attempts to control the roll angle in a normal manner 
also helped to reduce the excursions of the sideslip 
angle . 

A by-product of the computer study was the confirm- 
ation that the frequency- of the oscillation was not a 
significant contributor to the basic instability. Quite 
obviously, this does not mean that frequency can always 
be dismissed as insignificant, but only that an instability 
of the type studied could exist even at frequencies 
completely within the control limits of a human being. 

In this phase, it was found that motions with a period 
as high as five seconds were completely unmanageable 
as long as a normal pilot technique was used. 

The flight phase of this investigation yielded 
results which essentially complemented those obtained 
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from the fixed base simulator. This was in spite of the 
fact that the character of the oscillatory motion obtained 
with the NAvion variable stability airplane was consider- 
ably different from what was observed with the fixed 
base simulator. 

As outlined previously, the connection seemed to be 
that, in either case, an attempt by a pilot to control 
the roll angle by manipulation of the ailerons resulted 
in an unfavorable generation of sideslip. Inasmuch as 
the rolling and sideslipping motions were highly coupled, 
the increased sideslip generated additional rolling 
moments which in turn required even greater control 
movements. The result was a coupled^ rolling, sideslipping 
motion of ever increasing amplitude. That the oscillations 
could be contrdlled - or even damped - by judicious use of 
the rudder, substantiated the validity of the foregoing 
explanation of the instability mechanism. 

The result of the combined studies was that certain 
combinations of roll angle, sideslip angle, and side- 
slip resultant from an aileron deflection, were found 
to be unfavorable. In particular, it was noted that 
the pilot-induced unstable behavior occurred only when 
the pilot, in an effort to return to a wings level 
condition, could not reduce the bank angle without 
increasing the sideslip angle. Thus an "effective" 
proverse aileron yaw was found to have a destabilizing 
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effect if the Dutch roll oscillation v/as of such a nature 
that the bank angle and the sideslip angle were in phase. 
For an "effective" adverse yaw, the opposite phasing was 
unfavorable. Multiple exposure photographs depicting 
the two motions described above are included herein. 

(Figs. 18 and 19). The term "effective" as used refers 
to the sense of the initial sideslip produced by an 
aileron deflection. This initial sideslip depends on 
the combined effects of Cn&q , C n p and Oi Q , rather than the 
actual sign of the yawing moment produced. Thus it might 
be termed an integrated effect. 

In addition to the above stated phase relationships, 
certain conditions were necessary before the instability 
could occur. First, it was necessary that the lateral- 
directional motions be highly coupled. Physically this 
was indicated by the large oscillations which occurred in 
the time histories of roll rate response to a step aileron 
input* Secondly, it was necessary that. the Dutch roll 
mode be lightly damped. No formal proof of this last 
statement has been offered herein, but it seems quite 
obvious that a highly damped mode would virtually never 
demand pilot attention. Practically speaking, the pilot 
would probably never be aware of its existence. Accord- 
ingly, one could hardly expect the pilot to ever initiate 
the destabilizing mechanism described above. From a root 
locus standpoint, an increase in the Dutch roll damping 




Fig- l& 

Dutch Roll Motion 



J3 and of Opposite Phase 




* 






Fig. 19 

Dutch Roll Mo+i on 
/3 and (p of +he Same Phase 
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would shift the Dutch roll poles and zeroes (Pigs. 8a 



and 8b) to the left. This would obviously have a 
stabilizing effect. 

Prom the discussion of the destabilizing mechanism, 
it is apparent that the phasing between the bank angle 
( <p ) and "the sideslip angle ( ) is an important indicat- 
or of how a closed loop instability might develop. In 
particular, this phase relationship determines the sense 
of the effective aileron yaw that can be regarded as 
unfavorable. Accordingly, a study was undertaken in an 
effort to ascertain what parameters determined this phase 
relationship. This study is outlined in Appendix C. 
Briefly, the phase angle between <p and was found to 
be heavily dependent upon the sign of the term:- sycj/jg 
This same combination of derivatives has been found to be 



important in the determination of the 



CO«t 






parameter 



Further analysis led to the conclusion that the 
phasing of the Dutch roll in roll could be expected to 
bear some relationship to the character of the roll rate 
response. Figures 9 and 16 showed that a distinct 
difference existed between the characteristic time 
constants of the roll mode for the favorable and un- 
favorable situations. In particular, the time difference 
was seen to be approximately equal to one half of the 
Dutch roll period. This effect is discussed in detail 
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in Appendix D. Briefly, a shift of 180 degrees (Dutch 
roll in roll) was seen to coincide with a shift through 
unity of the ■ magnitude of the parameter. 

In summary, for the type of pilot induced instability 
outlined herein, the arguments presented appear to be a 
physical way of describing the phenomenon depicted by 
the root locus studies of Taylor. The theoretical 
effort conducted in conjunction with this investigation 
seems to show that the onset of an unfavorable combination 

of phase relationships is coincident with the parameter 

60 * 

described as becoming greater than unity. 

Although only a limited number of cases of pilot- 
induced instabilities were studied in detail in con- 
junction with this investigation, the authors have seen 
no reason to suspect that the destabilizing mechanism 
as outlined herein cannot be regarded in a broader sense. 
Thus it is the opinion of the authors that an explanation 
of this type might well be valid for other instabilities 
known to be a function of pilot gain. 



- 68 - 



SHCTIOK V 
OPK CLUSIOKS 

As a result of this investigation, the following 
conclusions were drawn: 

1. The NAvion variable stability aircraft can be 
used for a general study of pilot-induced 
lateral-directional instabilities . 

2. A pilot induced instability of the type investi- 
gated can be successfully demonstrated f and 
analyzed with a fixed base simulator incorpo- 
rating an electronic analog computer. 

3. For this type of pilot induced instability of 
the Dutch roll mode to occur, certain conditions 
were found to be required. These are listed in 
a qualitative fashion as: 

a. A light ly damped Dutch roll transient; 

b. Highly coupled rolling and sideslipping 
motions - as manifested by a large per- 
centage of Dutch roll in roll. 

4. If the conditions listed in item (3) were prevalent, 
certain phase relationships between the roll angle, 
sideslip angle , and the effective aileron y aw were 
found to result in an unstable closed-loop system. 

In particular, if the Dutch roll oscillation showed 
the roll angle and the sideslip angle to be of the 
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same sign, an unstable system resulted if the 
effective aileron yaw was proverse. For the 
opposite phasing, an effective adverse aileron 
yaw resulted in an unstable closed loop system. 

5. For the type of pilot induced instability 
studied, the frequency of the oscillation was 
found to be an unimportant parameter. 

6. Even a marginally stable system (open loop) 
could be managed quite easily if a direct means 
for controlling the sideslip angle were provided. 
Accordingly, a system which was closed loop un- 
stable with only aileron control was found to be 
manageable if the rudder was used Q Also, the jb 
technique could be used, although greater care 
was required. 

7. It was further concluded that this study provided 
a physical interpretation of the type of pilot- 
induced instability predictable by the method 
described by Taylor (Ref. 2). 
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TABLE I 

PHYSICAL CHARACTERISTICS OP THE NAvion 
WING DATA 



Total Area (including ailerons, flaps 




and 19.87 ft. covered by fuselage) 


183.34 ft. 2 


Span 


33.38 ft. 


Aspect Ratio 


6.05 


Taper Ratio 


0.54 


Dihedral Angle 


7.50 deg. 


Root Chord 


7.20 ft. 


Mean Aerodynamic Chord 


62.35 in. 


Incidence Angle 




Root 


2.0 deg. 


Tip 


-1.0 deg. 


Sweepback of Leading Edge 


3.0 deg. 


Twist 




Geometric 


3.0 deg. 


Airfoil Section 




Root 


NACA 4415R 


Tip 


NACA 6410R 


Flaps, 40 deg., Plain 




jERON DATA (For one aileron) 




Area 


2.16 ft. 2 


Span 


61.99 in. 



30 deg. up; 20 deg. dwn. 



Deflection 
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TABLE I (continued) 

AILERON DATA (continued) 

Control 

Aerodynamic Balance 

Static Balance (Outboard end of 

each aileron) 

Trim Tab (Right aileron) 

Ratio of Aileron Chord to 
Wing Chord 

HORIZONTAL TAIL DATA 

Total Area (Including 2.37 ft. 2 
covered by fuselage) 

Span 

Aspect Ratio 
MAC 

Airfcil Sections, Root and Tip 
Incidence Angle 

ELEVATOR DATA 
Total Area 
Span 

Deflection 30 deg. 

Deflection Trim Tabs 

(32 in. span, 4-1/2 in. Chord) 

Root Chord 
Tip Chord 



Wheel Throw 
Frise-type nose 

Streamlined weight 
Fixed Bend Tab 

0.284 

43.05 ft. 2 
13.17 ft. 

I 

4.02 
3.34 ft. 

NACA 0012 
-3.0 deg. 

14.10 ft. 2 
73.58 in. 
up; 20 deg. dwn. 

+30 deg. 

1.5 ft. 

1.0 ft. 
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TABLE I (continued) 

VERTICAL TAIL DATA 

Area (Including 2.57 ft. 2 blanketed by- 

fuselage and excluding 1.84 ft. 2 dorsal 

fin area) 12.93 ft. 2 

Span 4.05 ft. 

Airfoil Section 

Root NACA 0013.2 Modified 

Tip NACA 0012-64 Modified 

Incidence Angle (With respect 



to FRL) 

RUDDER DATA 
Area 

Deflection 
Trim Tab 

Rig (Angle with respect to 
fin center line) 



2 deg. Nose Left 

6.05 ft. 2 
17 deg. L; 23 deg. R. 

Fixed Bend Tab 

3 deg. right 



FUSELAGE DATA 

Fuselage Length Over-all 
Width, Maximum 
Depth, Maximum 
Fineness Ratio 



27.25 ft. 
4.14 ft. 
4.40 ft. 
6.2 



POWER PLANT DATA 

Airplane is powered by one Continental E-185 engine. 
Maximum continuous rated horsepower at sea level 
185 at 2300 RPM. 
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TABLE I (continued) 



POWER PLANT DATA (continued) 

The propeller is a Hartzel hydro-selective propeller 
with the following characteristics: 

Activity Factor 100 

Diameter 86 in. 

Pitch (.75R) 26.5 deg. 

MISCELLANEOUS DATA 



Weight 

Basic 

Fuel (40 gal.) 

Pilots (2 with parachutes) 

Gross Weight 

Center of Gravity Position 
Tail Length 

Autopilot Mechanical Gear Ratios 
. Deg. Elevator per Deg. Servo-drum 
Deg. Aileron per Deg. Servo-drum 
Deg. Rudder per Deg. Servo-drum 



2129 lbs. 
240 lbs . 
410 lbs. 
2779 lbs. 
29.5$ MAC 
16.88 ft. 

0.919 

0.910 

0.330 
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TABLE II 

SPECIFICATIONS OF INSTRUMENTATION COMPONENTS 



RATE GYROS 
Manufacturer 
Model Number 
Power Input 
Pov/er Consumption 

Rotor Speed 
Weight 

Maximum Turn Rate 

Potentiometer Resistance 
Potentiometer Excitation 



SIDESLIP VANE 
Manufacturer 
Model Number 
Potentiometer Resistance 
Potentiometer Active Angle 

Potentiometer Excitation 



Minneapolis -Honeywell 
JG 7005A-24 

115 volts, 400 cycle A.C. 

52 watts (starting), 

13 watts (running) 

20,000 RPM 
1.75 lbs. 

Modified to approximately 
30°/sec. 

530 ohms 

Gyro # 9 - 24 volts A.C. 
Gyros #10, #11 and #13 - 
30 volts A.C. 

Giannini 

2516 

1997 ohms 

Yaw channel - + 30° 

Roll channel- +45° 

Yaw channel - 15 volts A.C. 
Roll channel- 30 volts A.C. 



- 76 - 



TABLE II (continued) 



TELEMETER TRANSMITTER UNIT 
Manufacturer 
Model Number 
Input Signals 
Information Channels 

Sampling Rate 

Accuracy 

RF Power Output 

Frequency Range 

Frequency Stability 

Primary Power Requirement 

Weight 



AS COP 

DT-4 

0-5 volts D.C. 

43 plus two for synchro- 
nization of ground station 
20 RPS 

+ 1 # of full scale 
4 watts 
215-235 MC 
+ 0.05# 

28 volts D.C. 

19.6 lbs. 
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TABLE III 

Dimensional derivatives of the X-15 at Mach = 3.0, (X 0 = 10° 



L (3 


= 


+9.02 


N p 


-.00451 


L P 


= 


-.212 


N r = 


-.0451 


L r 


S 


+ .172 


N *a = 


+ .454 


n p 


= 


+2.75 


II 

O 


+8.40 






TABLE IV 






Non dimensional stability derivatives of 


the X-15 


at Mach = 


3. 


O 

9. 

0 

11 

0 

0 






% 


= 


+ .044 


Cn r = 


-1.40 


% 


= 


-.280 


II 

O 

c 

a 


+ .052 


% 


= 


+ .228 


C ^a = 


+ .041 


Cn j$ 


= 


+ .315 


C u = 


+ .440 


Crip 


= 


-.140 


Cy„ = 

t 


-1.38 



TABLE V 

Autopilot potentiometer feedback gain constants. 

k 4 = +.541 deg. rudder/deg. sideslip 

^5 = -.121 deg. rudder/deg. /sec roll rate 

kg = -.111 deg. rudder/deg. /sec yaw rate 

kg = +.609 deg. aileron/deg. sideslip 

kg = -.0756 deg. aileron/deg. /sec yaw rate 

k 7/^1 3 _ -.345 deg. rudder/deg. aileron 
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AI TBN-DIX A 

Calculations for the Plight Simulation : 

Non-dimensional stability derivatives chosen as 
typical for the X-15 at M = 3.0, 0 C o = 10 ° are: 



= + .044 


Cn p = - .140 


C N - C u = .44 


Cip = - .280 


C n r - - 1 . 40 


C , 6 - -1.38 


Cjj p = + .228 


C„ io = + .052 


C nfi = + .315 


Cn. = + .04-1 




For one - g flight: 

f ■ JU-- 

1/2 C L S 

13,445 







1/2 ( .44 )(3 x 1005 )^( 200 ) = 3.36 x 10 ~ 5 lb-sec 2 /ft 4 ( 1 ) 



Where 1005 ft/seo is the speed of sound corresponding to 
an altitude for which p = 3.36 x 10 ~^lb-sec 2 /ft 4 
q = 1/2 pV 2 = 1/2 (3.36x10" 5 )(9.1x10 6 ) = 153 lb/ft 2 

T= M/pSV = 3 ^ 3 6 xlO ~ 5 (200)(3015)(32.174) = 20,6 
pi.= M/pSb = 3^3^ 4 f 0 -5 (200) (32. 174) (22. 36) = 2790 
k 2 = lx = 3.348 x 32.174 = 8.02 ft 2 



M 13,445 



' Z ^445 



= 78,691 x 32.174 = 188.4 ft' 



- 2 1 ? 



k„ 2 



= 2 8.02 = .0321 

5M“ 



( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 



J z = 



= 2 188.4 = .753 

W” 



2 



kz 

,TT* / 



2 



(9) 
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Values for dimensional derivatives were obtained as 
follows i 



L(i = 


.044(2790) 


+ 9.02 


(10) 


1 T*J* 


424. 4( .0321) 






II 

* 

q O 

<=* i 

II 

a. 


-.280 


= -.212 


(11) 


2(20. 6)(. 0321) 




_ C* r _ 
r 2.TJ X ' 


.228 


= +.172 


(12) 


"2(20. 6) (.0321) 








2790C.315) = + 

T2T7tfTT53T 


2.75 


(13) 


Kip 


- .14 


= -.00451 


(14) 


p erJi 


5T2o.SK. 724) 


Nr ZTJ t ~ 


- 1.4 


= -.0451 


(15) 


2 (20.6)( .724) 

j 








.052(2790) 


+ .454 


(16) 


424.4 (.753) 






; - uCjl&a _ 


.041 (2790) 


+ 8.40 


(17) 


424.4 (.0321) 







The equations which express the relationship between the 
simulator airplane derivatives and the simulated deriva- 
tives are derived in Reference 4. The resultant express- 
ions for the autopilot feedback gain constants are: 



k 4 = 1.312 - 2.415 x 10“ 3 Lj 3 - 0.2774 Np (18) 
k 5 = -0.123 - 2.415 x 10" 3 ^ P - 0.2775 Np (19) 
kg = -0.124 - 2.415 x 10 -3 L r - 0.2775 N r (20) 
k Q = 0.346 + 0.0490 La, + 0.0278 No, (21) 
k g = 0.373 + 0.0490 L p + 0.0278 N p (22) 
k 1Q = 0.0828 + 0.0490 L r + 0.0278 N r (23) 



Values obtained, fron the foregoing equations, for feedback 
gain constants were: 

k^ = + .541 deg. rudder/deg. sideslip 

k^ = - .121 deg. rudder/deg/sec roll rate 

kg = - .111 deg. rudder/deg/sec yaw rate 

kg = + .609 deg. aileron/deg sideslip 

kg = + .562 deg. aileron/deg/sec roll rate 

k^Q = - .0756 deg. aileron/deg/sec yaw rate 
k^/k^ = - .345 deg. rudder/deg. aileron 

Cockpit potentiometer settings were obtained from the 
calibration curves for these values of feedback gain. 



- 81 - 



APPENDIX B 

The following is a derivation of the root locus 
equation used to study the effect of pilot gain. 

A pilot response,, proportional to roll angle, of 
the form K(1 + .57d) is assumed. 

The determinant set, in non-dimensional form, is 
as follows: 





P 






<$> 


da 






c 


v Zd 


-a 




C L + 2<x 0 d 


0 








° S 9 / Z 


?|Pd- Jxd 2- 


/i c i<s 


a 


(i) 


, u 


c Yz 


- J t d 


Qn£d 




>q 






0 


0 




K( 1 + .57d) 


-i 






Expansion of the above can he symbolized in determinental 


form as: 
















C^- 2d 




a 


Cl + Zaod 










yU Cj2|8 


Cq 


r/a 


5«£ _ J,d^ 
e. 










jmCnp 


Cnr - J £ d 

o *■ 


Cn £d 




i 


















(2) 








Cya - 2d ~Z 


o 








+ K(l + 


51 d) 


P<Cifb 


Ottr/g. 




= 0 










yuCnA 


Cnr /2 ■ J 4 d 


t-hCndo 
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Expansion of the first determinant set will yield 

a quart ic equation which can be divided by 2J J and then 

X z 

factored into a form from which the characteristic modes of 
transient motion can be deduced. Thus this is the character- 
istic equation and will be symbolized as: 

|D| = ( A + /Vr)(X + \ s )( \ ^ + 2.^<f>C0ni|j + OOnci) 2 ") , 



The remaining determinant can be expanded and the 
sum represented as follows: 

I D| * K(K57d) (mC| J t d) + ZfxCn^ 









uCn^Q [(CyjQ,- 2.d)^- + 



This can be rearranged to yield: 



I D | + K(l ^.57d)M Q<Sq jfa 2 -- (Cn r b J*C vft)d -t ^y^ + 2/qCnJ 

Jx ( 2 J*. 



(5) 



+ Cn&q. Gir r d / CyA + Zu Qtft'l 

C4i a 2JiL If" 'Cjr J 



The terms multiplying the human transfer function are 
seen to be a quadratic which is known to yield a complex 
pair of roots. Accordingly the above may be symbolized 
as follows: (d is replaced by A ) 

|D| + jutOiSq K(l + .57A)(V + + COn/) 

7*. 



( 6 ) 



The root locus equation is thus found to be: 

(A b Ar)(A+ As)( A 2- + 2 .^ IPncfA + OJm t) 2 '] + ]UtQl 6 q [7 = 0 
(\ z + Z .^<pCDn<pX + a) n 0 2 ')(l + .57A) Jx 
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APIENDIX C 

Discussion of Phasing of Boll and Sideslip : 

One of the parameters discussed in conjunction with 
the instability studied herein is the phase relation of 
roll angle and sideslip angle. This relationship is the 
subject of the following discussion. 

Prom the basic determinant for the lateral- 
directional motion of the airplane, given here in Laplace 
operator form as: 



I 5 

y uCn /5 



On. 



C# 

-Z 

CJr 

Z 

- 



<t> 



6a 



P 



yuCjlyb 



O 



.ym. Cn6a 



Cl + ZotoS 


0 




-yaCj?6a 




_ juGn<Sa 
' S 


(p to/3 as: 




-Z 


0 


Cfir 

C " r -J x 5 

Z. 


_ /lGn<Sa 


-Z 


C L + ZtXoS 


Ci r 

Z 

Cnr _ ] s 
z 


%S- JxS 1 
cw 



( 1 ) 



( 2 ) 



s 
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This may be written as: 



ft/n = S 2 + A s + B 

P 



C s + D s + E 

where 

- Gn6a C.£r G> 



A = 



^Y(2> _ Cnr 
Ci 6 q ZJz Z 



B - M^n a C^Cn n _ Cn8g Q(?r 't 1 y/3 , 

' Jt 4 j; Ciia £JZ vz Hjjfc 



G = Cn6q _ 2.0(o Ji 

Ci6a 



D _ Cnp - ClJ £ + CXoCrir — Gn8a (C£p + OcoCip) 

G(Z<ba 



(3) 



p _ GlCht _ Gncba CuCQr 

Z CYSo ~ 2T 



p 

For the condition of low dynamic pressure (153 lb/ft ), 
with CY 0 = 10° this ratio can be written 



<P 

P 



u = a 






+ 1.81s + 962.1 
"2 



or 



-.182s - .366s - .3716 



(4) 



s ± .91 - 31.0 .1 

s + 1.01 + .957 j 



Constant 



These roots are plotted on the complex plane and yield a 
phase relation for of -173 degrees (Pig. 20). 

It appears from an examination of the above expression 
of $ / jb that in those cases where the Dutch roll mode is 
lightly damped, the damping terms of the numerator and 



imaqinarv 

axis 




FIG. 20 

0/6 PHASE 
RELATIONSHIP 



real 

axis 



0//3 PHASE - Ih \ /A - ZL - Z2 

= -173° 



- 86 - 

denorainator will also be small. This means that all the 
roots will lie close to the imaginary axis. Further, 
by analyzing the plot of these roots, it is evident 
that the phasing of the 0/j3 motion is primarily de- 
pendent upon the relative position of the numei*ator 
root and the Dutch roll root. If the numerator root 
has an imaginary part greater than the imaginary part 
of the dutch roll root, the phasing between 0 and jb 
will be approximately 180 degrees. If the Dutch roll 
imaginary part is the larger, then the phasing will be 
nearly zero degrees. 

Those parameters which affect the imaginary part 
of the Dutch roll root are difficult to isolate since 
all approximations applied to the condition under investi- 
gation yielded poor results. However, it is evident 
that a major contribution comes from This 

combination of parameters also appears in coefficient B 
in the expression for 0/y3 , so its influence in shift- 
ing the two roots is equal. The only other significant 

term in coefficient B is -2.ju.C0a _Qn6q_ . it is through 

' 1 Cfi-ba 

this expression that a major shift in the numerator root 
is achieved. In the case under consideration, this one 
term is large enough to shift the numerator root below 
the Dutch roll root, changing the phasing to zero degrees. 
This shift would occur if one of the three parameters had 
a sign change. 
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This is merely additional supporting evidence for 
the argument that the phase relation between^ 5 and jb 
and the yaw due to aileron are important parameters 
in the explanation of pilot induced lateral-directional 
instabilities. It also points out the problem which 
will arise in any airplane - pilot combination where 
Cji/5 changes sign within the performance envelope. If 
this change of sign occurs in an area where other 
conditions are favorable to pilot-induced instabilities, 
then these instabilities may well be encountered. 
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APPENDIX D 



Discussion of Dutch Roll to Roll Mode Phasing : 

In an attempt to shed more light on the mechanism 
producing the type of instability under consideration, 
the following analysis of Dutch roll to roll mode was 
made. Of particular interest was the phase relation 
of the two motions. 

Again the basic determinant for the lateral- 
directional motion of the airplane is given here in 
Laplace operator form. 



p 


d<P 


<P 


6ct 


Cyj5“ 2 s 


-z 


Cl+ Zol 0 s 


0 


pCp 


CCr 

2. 


9§£s - JxS^ 


-pQ. 


uCnft 


J*s 


%’s 





Forming an expression for 0 in terms of this 
determinant and applying a step input yields the follow- 
ing determinant for fp . 

C.p-2s -2 O 

Or/Z -JxCjl&a 
Qn r _ J^s -/iC n 6a 



cp = 



Cyp 
u iCjija, 
^Cnjb 



-2 

CJlr/2. 



Z 



Cl 

- JxS 2 " 



Qnr-J £ S Qnp S 
2 a 



This may be written as: 



<t> 



a( 



Bs 



C) 



(s - X ,)(s - \zj[s - ~ X 4) 



- 89 - 



where 




q CruSg _ (Crv + 

2- Cj^q Z J t 

s' ( C yft C n n/Z.) -f 2-M.Cn/S _ CCr Cn 6o 

ZJ t ~ ZJ t C46a 

and the denominator represents the characteristic of 
motion. 

For the condition of low dynamic pressure (153 lb/ft^) 
with 0( 0 - 10°, this expression becomes: 

0 , 3565 (s 2 --)- I.8l4s 9b Q 

(s - X i j(s - Xi^s - Xa)(s - X^ 

The roots of the above are plotted on the complex 
plane (Fig. 21). It can be seen that the phase angle 
changes by 180 degrees whenever the numerator root moves 
inside the Dutch roll root. Further, a comparison of the 
roll traces (Figs. 9 and 16) shows the same resp.lt. 

Figure 9, with effective adverse yaw, shows a time delay 
of 1.5 seconds to reach the first peak on the roll trace. 
Figure 16, with the same conditions except effective 
proverse yaw, shows a time delay of .5 seconds. This 
difference in time corresponds to one half the period of 
the Dutch roll motion. Here again, as in the <fi to fd 
relation, we observe a change in phase of 180 degrees at 
the same point where 03^ /cOcjj becomes greater than 1. 
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